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INTRODUCTION 


A?i  objective  procedure  for  producinr)  a three-dimensional  mathematical  • 
model  of  sound  velocity  in  the  oceans  has  been  developed  by  tlie  U.  S. 

Naval  Oceano<iraphic  Office.  This  model  contains  the  necessaiy  statistical 
characteristics  to  form  the  basis  for  acoustic  province  delineation. 

The  mathematical  procedures  necessary  to  generate  a primary  three- 
dimensional  model  of  sound  velocity  containing  horizontal  spatial  frequency 
components  in  the  vertical  plane  are  Iteing  included  in  tiie  complilation 
of  such  products  as  NAVOCCANO's  ASW  Predicton  Area  Charts. 

The  accuracy  of  this  modeling  system  for  predicting  sound  velocity 
at  depth  betii/oon  ?00  m and  2450  ivi  at  any  latitude  and  longitude  viithin 
the  model  area  has  been  demonstrated  to  be  well  within  2 m/sec.  RMS.  The 
significance  of  tliis  modeling  concept  has  been  demonstrated  through 
rapid  computer  generation  of  acoustic  province  chai'ts  as  mentioned 
above,  as  well  as  a sound  velocity  envelope  for  each  delineated  acoustic 
province. 

The  validity  of  the  model  has  been  proven  in  the  Indian  Ocean 
(Arabian  Sea),  as  elsewhere,  where  good  agreement  was  found  between 
sound  velocity  profiles  generated  by  the  model  and  measured  cruise  data. 

The  mean  difference  between  the  model  and  the  measured  data  in  the 
Arabian  .Sea  was  0.43  rn/sec. 
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DESCRIPTION  OF  THE  SOUND  VELOCITY  MODEL 


An  objective  procedure  for  produciiuj  a tiireo-diinonsional  liiathetTiatical 
model  of  sound  velociti'  and  ott(er  oceanographic  parameters  has  been 
developed  by  the  U.  S.  Naval  Oceanographic  Office.  Tin’s  model  contains 
the  necessary  statistical  characteristics  to  fnnn  the  basis  for  the 
prediction  of  detection  ranges,  acoustic  survey  design,  and  acoustic 
province  delineation. 

The  model  contains  the  mathematical  procedures  necessary  to  generate 
a primary  three-dimensional  model  of  sound  velocity  containing  horizontal 
spatial  frequency  components  down  to  60  nautical  mile  wave  lengths  with 
400-meter  wave  length  resolution  in  the  vertical  plane.  Multiple  tests 
have  shown  this  model  accurate  within  2 inetcrs/second  RMS  in  predicting 
sound  velocity  at  depth  between  200  and  2460  meters  at  any  latitude  and 
longitude  within  the  model  area. 

This  model  functions  as  follows: 

1.  The  latitudes  and  longiLudes  of  the  corner  points  of  a polygon 
inscribing  the  area  of  interest  define  the  search  area  from  which  sound 
velocity  profiles  will  be  e'ttractod  from  the  On-Station  (OSTA)  data 
tapes.  These  tapes  are  a packed  listing  of  all  NODC  station  files 
through  1970  and  ,.oritain  sound  volcoity  coniputed  by  Wilson's  1960 
equation. 
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2. 


A 1°  mcrcator  grid  is  constructed  to  contain  the  search 


polygonal  area  and  provide  for  a 2°  overlap  into  adjacent  areas  which 
assures  cent  inunity  of  contours  into  adjacent  areas. 

3.  Each  OSTA  Station  is  tested  to  assure  conformance  to  input 
selection  criteria. 

4.  Accepted  stations  meeting  minimum  depth  requi ri^ments  but 
failing  to  reach  the  desired  maximum  saniple  depth  are  augmented  by  sound 
velocity  extrapolated  from  surrounding  data  to  the  desired  maximum 
depth. 

5.  Accepted  profiles  are  subjected  to  cubic  spline  interpolation 
to  yield  a final  profile  from  200  to  2450  meters  with  values  at  every  50 
meters  (46  levels)  at  eacti  1"  grid  point. 

6.  A seventh  degree  least  squares  orthogonal  polynomial  is  then 
fitted  to  each  final  profile  and  the  coefficients  retained. 

7.  The  resultant  sets  of  orthogonal  polynomial  coefficients  are 
fitted  into  a 30  x 30  minute  grid  by  employing  a three-phase  algorithm 
which  interpolates  r'andom  data.  Phase  one  partially  fills  the  desired 
grid  by  assigning  input  values  to  the  closest  grid  point.  If  more  than 
one  point  lines  near  a given  grid  point,  its  value  will  be  a weighted 
average  of  these  input  values.  Pliasc  two  uses  a modified  version  of  an 
alcjorithin  developed  by  Shepard  (1968)  to  expand  this  partial  grid  by 
completing  boundaries  and  adding  enough  interior  values  to  allow  splining. 
Phase  three  than  completes  the  g»  id  using  a bicubic-spline  algoritlim 

of  Davis  (1970)  applied  twice  and  averaged. 
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8.  The  gridded  orthogonal  polynomial  coefficients  througli  the 
fourtii  degree  are  prefitted  for  smoothing,  and  the  magnitude  of  the 


gradients  are  computed  and  contoured  to  yield  secondary  province  l.'Oundaries 
from  selected  values  of  the  gradients  of  the  coefficients. 

9.  Polynomial  curves  are  computed  from  the  resultant  linear 
correlation  coefficients  to  yield  minimum,  mean,  and  maximum  sound 
velocity.  The  geographical  [lositions  of  grid  points  and  the  related 
linear  correlation  coefficients  and  groups  of  the  latter  characterize 
the  average  and  type  shape  of  the  sound  velocity  profiles  for  each 
secondary  pi'ovince. 

10.  The  Navy  Interim  Surface  Ship  Sonar  Prediction  Model  (NISSM) 
computes  for  each  secondary  sound  velocity  province, the  probability 
of  detection  ranges  in  ki loydrds,and  minimum,  mean,  and  maximum  sound 
velocity  profiles  constructed  from  the  given  linear  correlation  coefficients 
as  wel 1 . 

11.  The  distribution  of  detection  range  variability  within  the  | 

secondary  provinces  is  analyzed  to  yield  the  final  selection  of  primary  | 

sound  velocity  provinces,  j 

MODEL  FOR  THE  AR,ABIAN  SEA  ; 

I 

An  example  of  the  output  of  this  model  is  demonstrated  in  Standard 
Navy  Ocean  Area  and  Region  IN-1  (the  Arabian  Sea)  where  one  can  delineate 
areas  with  similar  coefficients  and  therefore  with  similar  sound  velocity 
characteristics.  Figure  1 indicates  that  this  sea  can  be  divided  into 
six  primary  sound  velocity  provinces  including  the  Gulf  of  Oman  and 
Persian  Gulf.  Figures  2 througii  7 depict  minimum,  mean, and  maximum 

sound  velocity  structure  in  each  province.  should  be  noted  that  each 

province  has  a distinctly  different  stiucturc. 

j 
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A coiiipari son  of  sound  velocity  profiles  generated  by  the  model  for 


this  area  with  measured  cruise  data  indicates  good  agreement  (Coche, 

1974).  Sixteen  stations  occupied  in  tlie  Indian  Ocean  during  the  south- 
west monsoon  were  compared  to  the  model  and  the  mean  difference  was 
found  to  be  as  low  as  0.43  in/scc.  (01.4  ft/'sec.).  figure  8 (after  Cocke), 
which  illustrates  comparisions  of  cruise  data  wi tit  the  model,  demonstrates 
the  model 's  validity. 

Detailed  analyses  of  the  effects  of  the  northeast  and  southwest 
monsoons  on  the  sound  velocity  structure  in  tlie  Arabian  Sea  (and  Indian 
Ocean)  is  contained  in  Fenner  and  Bucca  (1972).  In  addition,  tabulations 
of  sound  velocity  at  100-metcr  intervals  below  1,000  meters,  as  well  as 
information  on  critical  depths  for  the  Arabian  Sea  are  available  in 
Audet  (1975). 

SUMMARY 

A mathematical  model  has  been  developed  which  can  predict  sound 
velocity  at  deptti  between  200  and  2450  meters  at  any  latitude  and  longitude 
within  the  model  area. 

Application  of  this  model  to  the  Arabian  Sea  resulted  in  the 
delineation  of  six  distinct  acoustic  provinces. 

Comparisions  of  sound  velocity  profiles  generated  by  tlie  model  for 
each  province  with  measured  cruise  data  indicate  that  the  mode)  is 
valid.  The  mean  difference  between  the  i’’OcIg1  and  measured  data  in  the 
Arabian  Sea  is  only  0.43  m/sec. 
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